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Abstract: This paper presents the formulation of a 
Lyapunov function for an exponentially stable linear time- 
varying (LTV) system using a well-defined PD-spectrum and 
the associated PD-eigenvectors. It provides a bridge between 
the first and second methods of Lyapunov for stability 
assessment, and will find significant applications in the analysis 
and control law design for LTV systems and linearizable 
nonlinear time-varying systems. 


I. INTRODUCTION 


Consider linear time-varying (LTV) system described by 

x = A(t)x, x(to) = Xq, t> to (1) 


where A(t) is continuous and bounded for all t > to . LTV 
system theory is an important tool for nonlinear system 
control analysis and synthesis. LTV systems arise not only 
from time-varying parameters and structures of the system, 
such as the changing fuel mass and staging of a launch 
vehicle, changing operating environment, such as the 
aerodynamic coefficients of an aircraft, they also arise from 
linearizing a nonlinear system along a time-varying 
trajectory which yields LTV error dynamics. 

LTV system control methods, such as linear matrix 
inequality control [1,2] and linear parameter varying (LPV) 
control [3] have been developed. A differential algebraic 
spectra theory (DAST) for LTV systems [4, 5, 6] extends the 
algebraic spectra concept of linear time-invariant (LTI) 
systems to LTV systems, which is a powerful tool to analyze 
and synthesize LTV systems. It is noted that the time- 
varying Parallel Differential (PD) eigenvalues defined in 
DAST is different from the frozen-time eigenvalues given as 
the roots of det[A(£) — A(t)] = 0. While they are related to, 
they are not the same as Lyapunov, Floquet or Bohl 
exponents. Compared to LPV and LMI, which is essentially 
a linear robust control approach, DAST grasps the intrinsic 
nature of time-varying dynamics of LTV systems, and 
potentially can result in less conservative controller design 
and improved performance. A nonlinear control approach, 
called trajectory linearization control (TLC), has been 
developed based on DAST and successfully applied to 
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challenging control problems [7-12, 23, 26, 27]. One unique 
feature of DAST based TLC is the ability to assign time- 
varying closed-loop dynamics for real-time adaptation of 
modeling errors, changing operating conditions and 
objectives, and unforeseen plant degradation and faults. 
Adaptive time-varying bandwidth (TVB) and time-varying 
notch filter (TVNF) schemes based on DAST have been 
developed to improve the TLC stability for actuator 
saturation and actuator failure, such as in the direct fault 
tolerant control developed in [11], as well as excitation of 
unknown or time-dependent structural modes [12]. 

DAST-based stability assessment belongs to Lyapunov 
first method. Lyapunov second method is more widely 
employed in control theory and applications. However, it is 
difficult to apply Lyapunov second method in LTV system 
analysis due to the difficulty in constructing a causal and 
bounded Lyapunov function for stable LTV systems. 

The LTV system (1) has a unique solution given an 
initial condition x(to ), which can be written as 

x(t) = $(t,t Q )x(to) (2) 

where 3>(£,£o) is called the state transition matrix (STM). It 
is well known [13, 14] that (1) is exponentially stable if and 
only if there exists a positive definite (p.d.) and decrescent 
quadratic Lyapunov function V (t , x) such that 

0 < ci||z ||2 < V(t, x) = x T P(t)x < C2 1 1 ^ 1 1 2 ? Vt 
where P(t) is bounded, symmetric p.d., i.e. 

P(*) = P T (*), 0 < c\I < P(t) <c 2 I 

and the time-derivative of V(t , x) is negative definite (n.d.) 
satisfying 

V(t,x) = —x T Q(t)x < — C 3 ||:e ||2 < 0 
where Q(t) is bounded, symmetric p.d., i.e. 

Q(t ) = Q T (t), 0 < C 3 I < Q(t) < c±I 
The positive constants Ck are given by 
Cl = nun (inf A fc (P(i))), c 2 = max (sup X k (P(t))^j (3) 
c 3 = nun(infA fc (Q(i))), c 4 = max(supA k(Q(t))j 

where A^(X(f)) denotes the kth algebraic (frozen-time) 
eigenvalue of the time-varying matrix X(t) as a root of the 
algebraic characteristic equation det[X(t) — X(t)I] = 0. 
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Knowing the bounds of P(t) is essential in LTV system 
analysis and synthesis. For example, the constants Ck can be 
used to estimate the envelope of the state trajectory by 

||a(t)|| < is||a(io)||e-*m (4) 


where k = 


and fi = which also asserts exponential 


stability by the definition. In linear adaptive control [15, 16] 
and nonlinear adaptive control [17, 18, 19], the (constant) 
matrices P and Q for the nominal exponentially stable LTI 
controller are used in the adaptive control law design to 
assure stability of the overall adaptive system. Those results 
can be extended to LTV plants if practical methods can be 
developed for constructing the time-varying P(t) and Q(t) 
matrices. 

There are two types of methods for finding time-varying 
P(t) and Q(t) matrices. The first method is by solving for a 
p.d. and bounded P(t) from the matrix differential 
Lyapunov equation (DLE). 


P(t) = - P(t)A(t ) - A T (t)P(t) - Q(t) (5) 

P(t 0 ) = Po = Pq, t>t 0 


for any given p.d. and bounded Q(t). It is well known that 
the DLE (5) has a unique bounded solution P(t) if the LTV 
system (1) is exponentially stable. However, this method is 
not practical as DLE for an exponentially stable linear (LTV 
or LTI) system is not stable forward in time, thus in general 
does not provide bounded p. d. solutions. One approach to 
construct P(t) is to solve the DLE using the reverse-time 
method by letting r = —t, rendering the DLE into an 
alternative form 


P(t) = P(t)A(t) + A t (t)P(t) + Q(t) (6) 
P(t 0 ) = Po = Pq, t>t 0 


If the solution reaches a steady-state, that steady state can 
then be used as initial values for the forward integration of 
(5). Reverse integration requires knowledge of A(t ) for all 
future time t > to, which is an anti-causal requirement. 
Thus, this method is effective only when A(t) has a 
predictable future, such as LTI system, linear periodic 
systems, finite-duration LTV systems with a A(t) that is 
predefined over the operating duration, and freezes to 
constant at the end of the mission. In particular, for a LTI 
system, DLE is reduced to an Algebraic Lyapunov Equation 
(ALE). 

PA + A t P + Q = 0 (7) 


In any case, this method is not applicable for adaptive LTV 
systems where the closed-loop dynamics are to be 
synthesized in real-time in response to unforeseen changes 
in the plant dynamics or operating conditions. The results of 
the Lyapunov function for both LTI system and LTV system 
are summarized in [20]. In [21], the solution bounds of a 
general DLE (6) are analyzed. In [22], a relatively 
conservative solution bounds for DLE (5) is analyzed. 

The second method constructs the P(t) for a given Q(t ) 
using the STM by [13, 14] 


/ oo 

$ T (r,i)Q(r)$(r,i)dr (8) 

Eq. (8) requires integration over the infinite future time- 
interval, which is anti-causal, and cannot be utilized in real- 
time control applications. Moreover, it requires the STM 
<fr(t, r) to be known for all future time, which is not 
practical to obtain for LTV systems even if A(t) is known 
for all future times, except for some special cases, such as 
the LTI systems. 

In summary, both approaches to construct P(t) for 
stable LTV systems are anti-causal and not applicable in 
real-time control application. This is not surprising as future 
behavior, such as exponential convergence as t — ► oo, of the 
state x(t) is determined by future values of A(t). The issue 
at hand is to find causally enforceable conditions, i.e. 
conditions that can be applied at the current time t on A(t) 
so that exponential stability will be guaranteed. This can be 
achieved using the D AST-based PD-eigenvalue assignment 
control [5, 6]. 

In this paper, Lyapunov function for exponentially 
stable LTV system is constructed using PD-eigenvalues and 
PD-eigenvectors, collectively called the PD-eigenstructure. 
Characteristics of P(t) and Q(t) are related to PD- 
eigenstructure of the LTV system. Most importantly, the 
constructed Lyapunov function is causal and applicable to 
real-time applications. 

The developed Lyapunov function is a useful tool not 
only as an alternative to the ALE (7) for solving P and Q in 
many Lyapunov control methods for time-invariant linear 
and nonlinear systems, such as adaptive control, but also for 
extending those methods to LTV and linearizable nonlinear 
time varying (NLTV) systems. Using the results developed 
in this paper, it is possible to assign the nominal controller 
as an LTV system, the causal P(t) and Q(t) matrix can be 
used in the adaptive law, as illustrated in [23]. Thus there is 
potential to design a less conservative and more robust 
adaptive control system. The main results of this paper also 
establish an alternative proof for the PD-eigenvalue stability 
criterion, along with estimate of the exponential envelope of 
the state trajectory. 

Following this introductory section, Section II presents 
a brief overview of the DAST to facilitate exposition of the 
main results, which are presented in Section III. Some 
illustrative examples are presented in Section IV. The paper 
is concluded with Section V. 

n. OVERVIEW OF THE DAST 

In this section we provide a brief overview of the DAST 
[4, 5, 6]. Consider the n-dimensional NLTV tracking error 
or operating error dynamics 

x = f(t,x ), x(t 0 ) = z 0 , t>t o (9) 

where / is continuous in t and differentiable with respect to 
x , and x(t) = 0 is the desired equilibrium state that needs to 
be stabilized. The time-dependence of / may be due to 
either time-varying parameters or time-varying nominal 
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tracking trajectory of the NLTV system. To begin with, we 
consider the special case where (9) can be linearized at the 
null equilibrium point (i.e. along the nominal trajectory) to a 
LTV system of the form 

x = A c (t)x , x(to) = xo, t> to (10) 

where A c (t) is continuous and bounded, and is in the 
companion canonical form 


0 



In— 1 

0 


—ot\ (i) 

-a 2 (t) ■■■ -a n (t) 


where J n _i is the (n — 1) x (n — 1) identity matrix. By 
letting y = x i, the state equation (10) is equivalent to an 
nth-order scalar LTV differential equation 

2/ (n) + oc n (t)i / (n_1) 4 1- oi 2 (t)y + oi\{t)y = 0 (11) 

y ( - k ~ 1 \t 0 ) = x 0k , k = l,...n t>t 0 


The general solution to LTV equation (11) can be written as 
y(t) = Y.C kVk {t) = C k eh n{T)dT (12) 

k = 1 k = 1 

where |i/ fc (i) = exp (j£/%(r) dr)}^ is a set of linearly 

independent solutions to (11), and {pk(t)}k=i * s ca ^ e d a 
parallel differential (PD) spectrum for (11) in which each 
pk(t) is called a PD-eigenvalue for A c (t). The PD- 
eigenvalues are solutions to an (n — l)th-order NLTV 
differential equation, called PD-characteristic equation 

A (p) = v;-\p) + a n (t)v;~ 2 (p) + (13) 

• * * + Oi^{t)'Dp(p) + Oi2(t)p + Oil {t) = 0 

where V p {p) =[f t + p(t)] ( p(t )) = p(t) 4- p 2 (t), and 
V k p = V p V k ~ x . In particular, for n = 2, the PD-character 
equation is given by the Ricatti equation 

p + p 2 + Oi2 {t)p H" ot\ (t) = 0 (14) 


Note that for real- valued coefficients ak(t) 9 the PD- 
eigenvalues pk ( t ) may be complex- valued in order to avoid 
singular solutions. Moreover, the synthesis formula for 
aift) from a PD-spectrum {pk(t)}k=i given by 


Oii(t) = pi(t)p 2 (t) + 


pi{t)p 2 it) - p 2 {t)p l {t) 


cx 2 (t) = —p\(t) — p 2 (t) — 


p 2 (t)-p 1 (t) 

P 2 {t)~ P\{t) 


(15) 


For each PD-eigenvalue Pi{t) 9 an associated PD- 
eigenvector Vi(t) for A c (t) is defined by 

A c (t)pi(t) - pi(t)pi(t) = p,i(t) (16) 


where pi(t) can be expressed using the PD-eigenvalue pi(t) 
as iM(t)=[ 1 V Pi {\} V%{\} ••• V n ~ l { 1}] T The 


diagonal matrix 

Y{t) = diag [p\ (t), p 2 (t),---, p„(t)] 

is called a PD-spectral canonical form for A c (t ), and the 
matrix consisting of all corresponding PD-eigenvectors 

M(t) = [/xi (i) | p 2 {t) | | p„(t)] 


1 

i 

i 

1) 

®*U) ■ 

• 0*(i) 

KW 

*£(!) ; 

• ®5L(i) 





is called a PD-modal matrix for A c (t ), and for the LTV 
system (10). A PD-spectrum (PD-spectral matrix) is said to 
be well-defined if the PD-eigenvalues pk(t) are bounded and 
continuously differentiable at least n times. A PD-modal 
matrix is said to be bounded diffeomorphic if it is 
continuously differentiable and invertible, with a uniform 
bound on M(£), M -1 (£), and M{t). A bounded diffeo- 
morphic PD-modal matrix constitutes a (stability preserving) 
Lyapunov transformation that relates the PD-spectral matrix 
T{t) and A c (t) by 

r(t) = M-\t)[A c (t)M(t) - M(t )] 

The above concepts can be extended to general LTV 
system (1) via a Lyapunov transformation that reduces A(t) 
to A c (t) 

A c {t) = L-\t)[A{t)L{t) - L{t)\ 

Then the PD-spectral matrix T(t) for A c (t) is called a PD- 
spectral matrix for A(t) 9 and the composite Lyapunov 
transformation V(t) = L(t)M(t) is called the 
corresponding modal matrix for A(t). A sufficiency 
condition for the existence of such a bounded 

diffeomorphism L(t) is that there exists a bounded, 
sufficiently differentiable n x 1 matrix B(t) such that 
{A(t), B(t)}is strongly controllable, i.e. 
inf)det[B(i) | V A B{t ) | ■■■ | V A n ~ l B{t ) ]| = c > 0, 
where V A B{t) = B(t ) - A(t)B(t), V\ = V A V k A ~ l 

The PD-eigenvalue concept is a natural extension of the 
eigenvalue concept for LTI systems in that, when the 
coefficients a k ( t ) = = constant, the PD-character 

(differential) equation (13) possesses constant solutions 
p k ( t ) = pk = constant that are conventional characteristic 
roots for (11), and eigenvalues for A c (t) = A c , for which 
(13) reduces to the conventional algebraic characteristic 
equation, and the synthesis formula (15) reduces to the 
conventional polynomial coefficient synthesis formula. 
Moreover, stability of the LTV system (10) or (11), and 
consequently the null equilibrium of the NLTV system (9), 
can be determined by the negativity of the real-part of the 
PD-eigenvalues. In particular, exponential stability is 
implied by R e(p&(£)) < — c < 0, for some constant c > 0, 
and for all k. 


3 



It is noted that, since {yk{t)}k=i are linearly 
independent, by definition PD-eigenvalues are always 
distinct. For instance, given a second order LTI system 


symmetric matrix, and r) is the state transition matrix 
for the LTV system z = A z (t)z. Let M(t) be a Lyapunov 
coordinate transformation and define 


x = 



1 

— 2A 


x , 


the system has repeated eigevalue A, whereas a PD-spectrum 
can be chosen as pft) = A and p 2 (i) = A + ^ for some 
constant d < to. It should be noted that PD-eigenvalues are 
functions of time. They are different from Lyapunov 
exponent, Floquet exponent or Bohl exponent, which have 
constant values. PD-eigenvalues are also different from the 
frozen-time eigenvalues of an LTV system. It is well known 
that frozen-time eigenvalues given as the roots of 
det[A(£) — A(t)] = 0 are in general neither sufficient nor 
necessary for stability assessment of a LTV system. Since 
evaluating a PD-spectrum for an arbitrary given matrix A(t) 
entails solving the (n — l)th-order NLTV PD-characteristic 
(differential) equation, which in general is difficult. Thus at 
the present time using PD-spectrum as an analysis tool 
remains a challenge, and requires significant future research. 
However, PD-spectrum as a design synthesis tool is not 
only practical, and very effective, as has been demonstrated 
in [7-12, 23, 26, 27]. The main results of this paper that 
follow bridge an important gap between the Lyapunov First 
Method based PD-spectral theory and the Lyapunov Second 
Method based control design techniques. 


x(t) = M(t)z(t) 

If the null equilibrium point for the LTV system z = A z (t)z 
is exponentially stable, then V x (t,x) = x*(t)P x (t)x(t) is a 
Lyapunov function for x = A x (t)x with V x (t,x) = 
—x*(t)Q x (t)x(t), where 

A x (t) = M~ l (t)A z (t)M{t) - 

P x (t) = M*-\t)P z {t)M~\t) 

Q x {t) = M*~\t)Q z {t)M-\t) ♦ 

Proof. Suppose the null equilibrium point for the LTV 
system z = A z (t)z is exponentially stable and let 

/ oo 

$*(T,t)Q z (T)$ z (T,t)dT 

It can be verified that V z (t, z) = z*(t)P z (t)z(t) is a 
Lyapunov function for z = A z (t)z with V z (t,z) = 
—z*(t)Q z (t)z(t) [13, pp. 155-159]. Let M(t) be a 
Lyapunov coordinate transformation and define 

x{t) = M(t)z(t) 


in. MAIN RESULTS 

The main result of this paper establishes a relationship 
between the first and second method of Lyapunov for 
stability assessment using the PD-spectral concept. 

Theorem 1. Let Y(t) = diag[pi(t), p 2 (^)r' Pn(t)\ be a 
well defined PD-spectral matrix for A(t) with a bounded 
diffeomorphic PD-modal matrix M ( t ) = [p\ ( t ) | /x 2 (i) | | 

p n {t)]> where pk{t) we the corresponding PD-eigenvectors 
for the PD-eigenvalues If 

Re{p k (t)) < -c < 0 (17) 


for some positive constants c, then the null equilibrium for 
x = A(t)x is exponentially stable. Moreover, V x (t,x) = 
x*(t)P x (t)x(t) is a Lyapunov Junction for x = A(t)x with 
V x (t,x) = —x* (t)Q x (t)x(t), where 

P x (t) = (18) 

and 

Q x (t) = —2 M* _1 (i) Re(Y(t)) M -1 (t) (19) 

where * denotes the conjugate transpose. ♦ 

Proof of Theorem 1 relies on the following lemma, 
which is an important result in its own right. 

Lemma 1. Let 


/ oo 

0* z (T,t)Q z {T)0 z (r,t)dT 


where Q z (t) is a continuous, bounded, positive definite, 


Then the null equilibrium for the LTV system x = A x (t)x 
is also exponentially stable, where 

A x (t) = M -1 (t)A z (t)M (t) - M~\t)M(t) 

Let 

P x (t) = M*-\t)P z (t)M~\t) 

and 

Qx(t) = M*~\t)Q z (t)M-\t) 

Note that 

z(t) = M~ l [t)x{t) 


Then 

P x (t) = M-\t)P z {t)M~\t) 


/ OO 

^ z {T,t)Q z {T)^ z {T,t)dTM-\t) 

/ OO 

/ OO 

&*x(T,t)Qx(T)$x(T,t)dT 


It follows from [13, pp. 155-159] that V x (t,x) = 
x*(t)P x (t)x(t) is a Lyapunov function for x = A x (t)x with 
v x (t, X ) = -x*(t)Q x (t)x(t). □ 

Proof of Theorem 1. Suppose the PD-spectral matrix T(t) 
satisfies (17) and the PD-modal matrix M(t) is bounded and 
diffeomorphic. Let z(t) = M~ l {t)x{t). Then z = Y(t)z. 
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Let Q z (t ) = -2R e(T(£)) > 2 cl > 0. It is readily verified 
that P z = I satisfies the Lyapunov equation 

p z = -P z T(t) - r*(t)p* - Q z (t) 

Thus, V z (t : z) = z*(t)P z (t)z(t) = z*(t)z(t) is a Lyapunov 
function for z = T(t)z with V z (t,z) = 2z*(t)RG(Y(t))z(t) 
< 0. By Theorems 4.10 and 4.12 of [13], the null 
equilibrium of z = Y(t)z is exponentially stable. Since 
M ( t ) constitutes a Lyapunov transformation, the null 
equilibrium for x = A(t)x is also exponentially stable. It 
then follows from Lemma 1 that V x (t,x) = x*(t)P x (t)x(t) 
is a Lyapunov function for x = A(t)x with V x (t, x) = 
—x*(t)Q x (t)x(t), where P x (t) and Q x (t) are given by (18) 
and (19), respectively. □ 

Remark 1. The proof of Theorem 1 provides an insight into 
the relationship between Lyapunov first method and second 
method for LTV systems. In particular, in the PD-spectral 
coordinates z = Y(t)z , where the modes of the LTV 
systems are decoupled, the Lyapunov function V z (t,z) = 
z*(t)P z (t)z(t) = z*(t)z(t) represents the total energy in the 
modes, and V z (t, z) = 2z*(t)RG(Y(t))z(t) gives the energy 
decay rate in terms of the PD-spectrum, which is twice the 
exponential decay rate of the envelope of the trajectories. 
The energy level and decay rate are “scaled” by the 
condition number of the PD-modal matrix in the original 
coordinates for the LTV system x = A(t)x. 

Remark 2. For LTI systems with distinct eigenvalues, the 
conventional (algebraic) eigenvalues and eigenvectors 
define a constant PD-spectrum and a constant PD-modal 


Remark 3. The technical difficulty for obtaining a causal 
expression of a Lyapunov function for an exponentially 
stable LTV system is circumvented by requiring the real- 
part of PD-spectrum to be bounded by a negative number. 
This result allows time-varying closed-loop dynamics to be 
synthesized for adaptation to unforeseen modeling errors 
and adverse operating conditions as they arise, thereby 
achieving best tradeoff between the closed-loop system 
performance, robustness and control energy consumption. 
Remark 4. For LTI systems, Theorem 1 provides an 
alternative method to the ALE for computing a Lyapunov 
function. Different P matrix can be obtained using different 
modal matrices, which amounts to choosing a different 
coordinate for the state space. 

Remark 5. Theorem 1 also provides an alternative proof for 
the PD-eigenvalue stability criterion. 

IV. ILLUSTRATIVE EXAMPLES 

In this section we use second-order LTI and LTV system 
examples to illustrate the utility of the main results, and to 
gain some insight into the relationship between the first and 
second methods of Lyapunov. We will also demonstrate the 
fallacy of frozen-time eigenvalue concept for LTV systems. 

4.1 LTI Examples 

In this subsection we consider some LTI examples with 
distinct eigenvalues, which will be denoted A&. 

Example 1 (The Companion Form, Distinct Real 
Eigenvalues). Suppose A is in the companion canonical 
form 


V(x) = x*(t)Px(t) and on 
V(x) = —x*(t)Qx(t) given by (3) can be used to estimate 
the envelope of the state trajectory by (4). Notice that, by 
letting P = M* -1 M -1 , the eigenvalues A k{P) = o|(M -1 ), 
which are the singular value squared for M -1 . Thus, the 
estimate of the peaking bound k = \J c^lc\ is actually the 
condition number of the (PD) modal matrix M, which is an 
indication of how singular the modal matrix is. In particular, 
the condition number attains its minimum when the modal 
matrix is unitary, in which case the peaking bound is k = 1. 
Thus, the peaking phenomenon is determined by the 
orthoganality of the PD-modal matrix, rather than by the 
exponential decay rate. Moreover, when M is unitary, 
M" 1 = M\ then P = = 7, and A k (Q) = 

— 2Afc(R e(T)). Thus, the exponential decay rate 

[i = — C3/(2c 2) = — minfc|Re(Afc(T)|, which is decay rate of 
the slowest mode. It follows that when two eigenvectors are 
nearly collinear, not only peaking of the trajectory will be 
large, the decay rate will also be slower than the single 
slowest mode. In practice, if there are design freedoms to 
assign the eigenvectors, they should be assigned as 
orthogonal as possible. These observations also apply to 
LTV case with a suitably defined induced matrix norm, such 
as || j 4( • )|| 2 = supt> io || ^4 (t) || 2 . 


Lyapunov 

function 

A = 

0 1 


0 

1 

the 

derivative 


—a\ —a ,2 


— A1A2 

Ai + A2 


Then the modal matrix is given by 
M 


1 1 
Ai A2 


The P, Q matrices are given in terms of the eigenvalues by 




A?+Ai 

— (A1+A2) 

p = 


(a 2 -ai ) 2 

(A2-A1) 2 



— (A1+A2) 

2 



(A2— Ai ) 2 

(Az-Aj ) 2 J 



2AiA2(Ai+A2) 

4 A 1 A 2 

Q = 


(A2 — Ai ) 2 

(A2-AO 2 


4 A 1 A 2 

— 2(Ai4-A2) 



(A2-A1) 2 

(Az-Aj ) 2 


Note that in this case, M cannot be orthogonal when both 
eigenvalues are of the same sign. Thus, peaking is 
unavoidable. For a\ = 2, a^ = 3, we have Ai = — 1, 
A 2 = -2, and 


P = 


5 3 
3 2 


, Q = 


12 8 
8 6 


with the singular values cr(M -1 ) = {2.6180, 0.3820}, 
condition number k(M) = 6.8541, and the eigenvalues 
A(Q) = {17.5440, 0.4560}. Thus the estimated trajectory 
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envelope is 

||s(i)|| < 6.8541 ||ar(i 0 )l|e _0 - 0333(i-to) 

which should be compared to the slowest mode e l<J K # 

Example 2 (Orthonormal Modal Matrix, Distinct Real 
Eigenvalues). Consider 


A = 


-1.5 0.5 

0.5 -1.5 


with eigenvalues Ai = —1, A 2 = —2. The modal matrix is 
given by 


M = 


0.7071 0.7071 

-0.7071 0.7071 


Note that M is orthonormal, and 

p= ; : , q = 


12 8 
8 6 


with = {1, 1}, k(M) = 1, and A(Q) = {4, 2}. 

Thus the estimated trajectory envelope is 

which should be compared to Example 1 above. * 

Example 3 (The Companion Form, Complex Eigen- 
values). Suppose A is in the companion canonical form 


A = 


0 

1 


_-ai 

-a 2 _ 



0 


Then the modal matrix is given by 

1 1 
<J + jus a — juS 


M = 


The P, Q matrices are given in terms of the eigenvalues by 

-<r(g 2 -|-fa> 2 ) 



■ (<r 2 +^) 

—a 


P = 

2W 2 

—O' 

2J 2 

1 

,Q = 


. 2J 2 

2J 2 . 

L 


o 2 


w 2 

—o 

J 2 J 


= -2 aP 


2 

1.4142" 

, Q = 

"2.8284 

2 

1.4142 

2 

2 

2.8284 


||z(£)|| < 2.4142 ||:c(to)||e 


-0.1213(*-t 0 ) 


A = 


-0.7071 0.7071 

-0.7071 -0.7071 


with eigenvalues Ai )2 = -0.7071 ± jO. 7071. The modal 
matrix is given by 


M = 


which is unitary, and 


0.7071 0.7071 

0.7071j — 0.7071 j 


P = 


1 0 
0 1 


, Q = 


1.414 0 

0 1.414 


with a(M *) = {1, 1}, k(M) = 1, and A(Q) = 
(1.414, 1.414}. Thus the estimated trajectory envelope is 

\\x(t)\\ < ||a:(to)||e- 7071 ^ 

which should be compared to Example 3 above. * 

4.2 LTV Examples 

In this subsection we first consider a LTV example that was 
developed to demonstrate the fallacy of frozen-time 
eigenvalue concept for LTV systems defined by the roots of 
det[A(i) — X(t)I] = 0. We will compare the Lyapunov 
function derived from the PD-eigenvalue concept and the 
fallacious frozen-time eigenvalue concept. A third LTV 
example is given to illustrate application of the main result 
to a time-varying bandwidth (TVB) filter. 

Example 5 (Parameterized Marcus-Yamabe Equation). 

Consider the parameterized Marcus-Yamabe (M-Y) 
Equation (1) with [24] 


A(t) = 


a — + acos(urt) \uj — asin(utf) 

— — asin(utf) a — — acos(utf) 


which was developed to show the failure of the sufficiency 
of the stability assessment using the negative real-part 
criterion using the frozen-time eigenvalue. The PD-spectral 
and PD-modal matrices are given by 


Note that in this case, M cannot be unitary. Thus, peaking is 
unavoidable. For a\ = 1, a 2 = \/2, we have 
A 1)2 = -0.7071 ± j0.7071, and 


T(t) = diag[2a 2 
M(t) = 


us. 


-\uj\ and 


cos {^ust) sin(^o^) 
— sin(^a^) cos(^utf) 


with cr(M -1 ) = (0.7654, 1.8478}, k(M) = 2.4142 and 
A(Q) = {0.8284,4.8284}. Thus the estimated trajectory 
envelope is 


Note that the PD-eigenvalues are constant and real, and 
exponential stability will be achieved for us > 4a > 0. Note 
also that the PD-modal matrix is orthonormal for all t. The 
P(t) and Q(t ) matrices are given in terms of the parameters 
a, us by P(t) = I and 


Q(t) = 


-2 a + us — 2acos(ust) 
—2asm(ust) 


—2asin(ust) 

-2a + us + 2acos(ust) 


which should be compared to the envelope defined by the 
real-part of the eigenvalues e -0 - 7071 (* -to ). * 

Example 4 (Unitary Modal Matrix, Complex Eigen- 
values). Consider 


The eigenvalues for Q(t) are given by { us , —4 a us}. Thus 
Q(t) > 0 for us > 4a, in agreement with the PD-spectral 
stability condition. The parameterized M-Y equation 
specialized to the original M-Y equation [13, Example 4.22, 
p.158] when us— 2 and a = 0.75, for which the PD- 
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eigenvalues are p\(t) = 0.5, p 2 {t) = — 1. By the PD-spectral 
stability criterion, the system is unstable. This is also 
verified with the fact that for these values Q(t) < 0 for all i, 
andV(t, x) = —x*(t)Q(t)x(t) > 0, for all t. Thus, any 
initial value x(to) will move away from the null equilibrium 
x = 0. Now, consider the frozen-time eigenvalues for the 

original M-Y equation, which are Ai } 2 (t) = —\±j^ 9 
which by the fallacious frozen-time eigenvalue stability 
criterion would imply exponential stability. The frozen-time 
modal matrix can be easily found to be 


Afrr(t) = 


1 — f sin(2i) 
-fcos(2i) + i^ 


1 — f sin(2t) 

— §cos(2f)-*^ 


This would lead to a positive definite Q(t) > 0 with 

eigenvalues \{Q{t)) = { 4 - 3 s1n(2t) . 7 [ 4 - 3 sin( 2 *) ] }> which 
would also (fallaciously) imply exponential stability. 
However, with the correct time-varying coordinate trans- 
formation and with the frozen-time eigenvalues treated as 
PD-eigenvalues, this stable system would actually be given 
by x = B(t)x , where 


B(t) = (MFr(i)diag[Ai(t), A 2 (t)] + M FT (t))M Fr 1 (f) 


8-6sm(2t)+9sm(4£)+24cos(2£) 

24sin(2i)-32 

|sin(2t) - I 


|sin(2i) + 1 

23-48sm(2t)+9cos(4t) 

24sin(2i)-32 


A 


Example 6 (A time-varying bandwidth filter). Consider 
an LTV system (1) with 


A(t) = 



2 (u n (t) 


where 0 < ^ < 1 is damping ratio and uj n (t) > 0 is the time- 
varying bandwidth. The PD-spectral and PD-modal matrices 
are given by T(t) = diag[pi (t), p 2 (i)] and 


M(i) = 


1 

Pi(t) 


1 

P2{t) 


where p h2 {t) = -(u n (t) ± jy/l - C 2 ^n(i)- 


m 2 ( 1 -c 2 K (<) 

Q(t) = 2 (u n (t)P(t) 


Wn(t) 

c 


c 

1 

u n (t) 


with 


A , U(t) + 1) ± ^ 2 (i)-l) 2 + 4C 2 ^(i) 

( ( )) 4c4(t)(l-C 2 ) 


K,(M(t)) = sup 


i>io2iU n (i)v/l - C 2 




and the exponent of the norm-bound on the trajectory 
envelope is given by 


_£3_ _ inf 

M 2c 2 i>io K 2 (M(i)) 


Figures 1-3 illustrate a numerical example with £ = 0.707 
and u n (t) as shown in Figure 1. Figure 2 and 3 show the 
elements of P(t) and Q(t). Figure 4 contains a plot of 
K(w n ) an d a Pl°t Of * 


V. CONCLUSION AND FUTURE WORK 

In this paper a method has been developed for constructing a 
quadratic Lyapunov function for an exponentially stable 
LTV system using a PD-eigenstructure for the system. The 
result establishes a relationship between Lyapunov first and 
second methods, providing explicit upper and lower bounds 
on the Lyapunov function using the singular values of the 
inverse PD-modal matrix, and an upper bound on the time- 
derivative of the Lyapunov function using the (negative) 
real-part of the PD-eigenvalue having the smallest 
magnitude (the slowest mode). An estimate of the 
exponential envelope of state trajectory is given in terms of 
condition number of the PD-modal matrix, which 
determines the peaking, and using the real-part of slowest 
PD-eigenvalue and the largest singular value of the PD- 
modal matrix, which determines the exponential decay rate. 
Illustrative examples are given to show the utility of the 
result, and to demonstrate fallacy of the frozen-time 
eigenvalue based stability analysis for LTV systems. 

The main result of this paper can be used in further 
development of the differential algebraic spectral theory 
(DAST) and DAST-based control techniques such as the 
trajectory linearization control. It can also be used to extend 
many known Lyapunov control techniques for time- 
invariant systems to time-varying systems. Moreover, it can 
be used to develop real-time implementable theoretic 
stability metrics for adaptive control of LTV and 
linearizable nonlinear time- varying systems. 
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